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Abstract: In supersymmetric standard models with multi Higgs doublet fields, selfcou- 
pling constants in the Higgs potential come only from the D-terms at the tree level. We 
investigate the decoupling property of additional two heavier Higgs doublet fields in the su- 
persymmetric standard model with four Higgs doublets. In particular, we study how they 
can modify the predictions on the quantities well predicted in the minimal supersymmetric 
standard model (MSSM), when the extra doublet fields are rather heavy to be measured 
at collider experiments. The B-term mixing between these extra heavy Higgs bosons and 
the relatively light MSSM-like Higgs bosons can significantly change the predictions in the 
MSSM such as on the masses of MSSM-like Higgs bosons as well as the mixing angle for 
the two light CP-even scalar states. We first give formulae for deviations in the observables 
of the MSSM in the decoupling region for the extra two doublet fields. We then examine 
possible deviations in the Higgs sector numerically, and discuss their phenomenological 
implications. 
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1 Introduction 

Although the standard model (SM) has enjoyed a marvelous success in explaining phe- 
nomena observed at collider experiments, the physics of electroweak symmetry breaking 
remains unknown. The experimental detection of the Higgs boson is the most important 
issue to confirm the standard picture for the origin of mass of particles, and the Higgs boson 
search is underway at the Fermilab Tevatron and the CERN LHC. On the other hand, the 
Higgs sector of the SM is known to cause the hierarchy problem [1] due to the quadratic 
divergence in the one-loop correction to the mass of the Higgs boson. In addition, there 
are several phenomena confirmed by the experiments which cannot be understood within 
the SM, such as the neutrino oscillation, the existence of dark matter and the baryon 
asymmetry of the Universe. Therefore, the SM must be extended to solve these problems. 

Supersymmetry (SUSY) is expected to be a good candidate of new physics. It can 
solve the hierarchy problem by the consequence of the nonrenormalization theorem [2]. 
The stabilized Higgs boson mass makes it possible to directly connect the electroweak 
scale with very high scales such as the Planck scale or that of grand unification. Super- 
symmetric extensions of the SM with the R parity also provide dark matter candidates [3] . 
In addition, various mechanisms of generating tiny neutrino masses [4-8] as well as those 
of baryogenesis [9-11] may also be compatible to supersymmetric models. 

The minimal supersymmetric standard model (MSSM) is a SUSY extension of the 
SM with the minimal number of particle content. It requires two Higgs doublet fields for 
cancellation of anomaly. The most striking phenomenological prediction of the model is 
that on the mass (m/i) of the lightest CP-even Higgs boson h. It can be calculated to be 
less than the mass of the Z boson at the tree level. Such an upper bound on rrih comes 
from the fact that the interaction terms in the Higgs potential are given only by D-term 
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contributions which are determined by the gauge coupUng constants. At the one-loop level 
the trilinear top- Yukawa term in the superpotential gives a significant F-term contribution 
to nifi [12-14], by which m/j can be above the lower bound from the direct search results 
at the CERN LEP experiment [15]. The calculation has been improved with higher order 
corrections [16-18]. Apart from nih, the masses of H, and the mixing angle a are a 
function of only two input parameters at the tree level; i.e., uia and tan/3, where rriA is 
the mass of CP-odd Higgs boson H is the heavier CP-even Higgs boson, are the 
charged Higgs bosons, tan j3 is the ratio of vacuum expectation values (VEVs) of the two 
Higgs bosons and a is the angle which diagonalizes the CP-even scalar states. In particular, 
there is a simple tree- level relationship among the masses of H^, A and the W boson 
as w-^^ = m\ + m^, where 'mu± and miy are respectively the masses of and W^. 
These characteristic predictions can be used to confirm the MSSM. 

If by experiments the Higgs boson is found and its mass turns out to be higher than 
the upper bound predicted in the MSSM, we must abandon the model. It however does not 
necessarily imply that the SUSY itself is excluded. For example, in the next-to-minimal 
supersymmetric standard model (NMSSM), where a gauge singlet chiral superfield S is 
added to the MSSM, the upper bound on mh is greater than that in the MSSM due to the 
F-term contribution from the trilinear XH1H2S interaction in the superpotential [19, 20], 
where Hi and H2 are two chiral superfields for the Higgs doublets. The effect of this 
term can gain the mass upper bound, so that nih can be as large as 140 GeV (400 GeV) 
by assuming that the running coupling constant A does not blow up below the grand 
unification scale ~ 10^^ GeV (the TeV scale) around tan/3 ~ 2 [20]. The similar effect on 
mil also appears in the model with additional triplet chiral superfields with the hypercharge 
of y = or y = ±1 [21, 22]. The other possibility to change the MSSM prediction on 
m/i is to introduce extra gauge symmetries which are broken spontaneously at the TeV 
scale. The mixing of gauge bosons of these new symmetries with those of the SM gauge 
symmetries yields the new D-term contribution to the Higgs potential, which can change 
for instance ruh [23, 24]. 

In general, interaction terms of the Higgs potential in extended SUSY standard models 
are composed of the D-term and the F-term as well as the trilinear soft-breaking term, so 
that extended SUSY Higgs sectors predict different masses and coupling constants from 
those in the MSSM. If the extra scalar fields are sufficiently light, they largely mix with 
the MSSM-like Higgs bosons, and consequently predictions on the MSSM observables are 
modified from the MSSM values at the tree level. In addition, such light extra scalars may 
directly be detected at the experiments. In this case the model can be easily distinguished 
from the MSSM. On the other hand, extra scalar bosons may be too heavy to be measured 
directly. Even in such a case the effect of these extra fields can change predictions on masses 
and coupling constants for the MSSM-like scalar bosons significantly at the low energy. 
This is not contradict with the decoupling theorem by Appelquist and Carazzone [25]. In 
fact, in the large mass limit of all the fields other than those in the SM, the model can be 
continuously reduced to the SM. 

In this paper, we consider the model (4HDSSM) where two extra doublet superfields 
are added into the MSSM in order to examine the decoupling property of extended Higgs 
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sectors without interactions from the F-term at the tree level. In this model, if there are 
mixings between the MSSM-like doublet fields and the extra doublet fields due to a large 
soft-breaking B-term mixing, nonvanishing effects can appear in the MSSM observables 
such as nih, mn, and a. Notice that this nonvanishing effect is not the so-called 

nondecoupling effect. The usual nondecoupling effect, which for example appears in the 
radiative corrections to the gauge boson two-point functions, does not vanish in the SM-like 
limit with taking the extra particles to be heavy. On the other hand, the present effect 
due to the large soft-breaking B-term mixing on MSSM observables in the 4HDSSM can 
be significant in the MSSM-like limit but does decouple in the SM-like limit by taking 
ruA — ?• oo according to the decoupling theorem [25]. Therefore, we call such a nonvanishing 
effect on the MSSM observables as the quasi-nondecoupling effect. We deduce analytic 
expressions of such quasi-nondecoupling effects on the MSSM observables when the extra 
doublet fields are heavy. Modifications from the MSSM predictions by these effects are 
then studied numerically. We find that the quasi-nondecoupling effect due to the B-term 
mixing can be significant. For example, ruh can be larger than about 10 % as compared to 
that of the MSSM prediction with the same MSSM input parameters, while m/^± can be 
about 20 % (10 %) smaller than the MSSM prediction when niA = 150 GeV (200 GeV). 
The mixing angle sin(/3 — a) also receives a significant deviation from the MSSM value. The 
Higgs potential and the mass matrices in the 4HDSSM have been analyzed by Gupta and 
Wells [26], but they have not included these quasi-nondecoupling effects discussed here. 

Suppose that the candidates of SUSY partner particles and extra Higgs bosons in the 
MSSM (ff, A and H^) are found at the LHC or the International Linear Collider (ILC) 
in future and that their properties look like consistent with the MSSM. In such a case, 
to confirm whether it is really of the MSSM or not, the MSSM relations are tested as 
accurately as possible by experiments. Our study can be particularly important in this 
case. 

In general, in extended Higgs sectors with multi-doublets, there necessarily appears 
the flavor changing neutral current (FCNC), which is strictly constrained by experiments. 
A simple prescription to avoid it would be imposing the (softly-broken) discrete Z2 sym- 
metry to the model [27]. We classify the type of Yukawa interaction under the discrete 
symmetry according to the assignment of the Z2 charges. Phenomenology of the 4HDSSM 
has been discussed by Marshall and Sher in the case with the lepton specific Yukawa in- 
teraction [28]. We here do not discuss influence of the heavy extra two doublet fields on 
the flavor physics, just assuming that the tree- level FCNC is forbidden by a Z2 symmetry. 
We then concentrate on the deviations in the Higgs sector due to the quasi-nondecoupling 
effects of these extra doublets. 

This paper is organized as follows. In Sec. II, the 4HDSSM is defined and the mass 
matrices for the Higgs bosons are obtained. In Sec. HI, we give a short description of 
the decoupling property in extended SUSY Higgs sectors, where the (quasi-)nondecoupling 
effect is discussed. We then present formulae for the quasi-nondecoupling effects on the 
MSSM observables. Numerical studies are also shown, and the results are discussed. Con- 
clusions are given in Sec. IV. We show the method of obtaining the Higgs potential in the 
useful basis from the general basis in Appendix A. The flavor structure in the 4HDSSM is 
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shortly discussed in Appendix B. 



2 The model 

We here discuss the 4HDSSM, in which two extra isospin-doublet chiral superfields H^, 
{Y = -1/2) and {Y = 1/2) are introduced to the MSSM in addition to the Higgs 
doublets Hi (y = —1/2) and H2 (Y = 1/2). The general expression for the superpotential 
with the R parity is given in terms of chiral superfields as 

W = - {Yu\jU'k,H2 ■ Qlj + {YdhD'jiiHi • Qlj + {Ye)ijE%^iHi ■ Llj 

- {YLhUkH, ■ Qlj + {Y^hD'j,,Hs ■ Qlj + {Y:hE'k,H, • Llj 

— Ml2-ffl • H2 — • H4 — ^32i?3 • H2 — fJ^S^H^ • J (2.1) 

where QLi, U'^^ (-^ra) are chiral superfields for the i-th generation left-handed quark doublet 
and right-handed up- type (down- type) quark singlet while Lu and E'^- are those for the i-th 
generation left-handed lepton doublet and right-handed charged lepton singlet, respectively. 
The most general holomorphic soft-SUSY-breaking terms with the R parity is 

^soft =-\iMiBB + M2WW + M3GG) 

{Ml\M^l + (M2)^34,t$3 + {Ml)\^^^, + (M2)33|.t$3 
{Ml)22^l^2 + {MlU^l^^ + (M|)24^'i^2 + (M|)44^I<^4 

+ {A,re*RA ■ hj + {Kre*R^^3 ■ + h.c. 

- [Bi2m^l ■ ^2 + B34m^3 ■ ^4 + BmUm^i ■ 'I4 + 832^^3 ' ^2 + ^.C.^ 

(2.2) 

where B, W and G are gauginos corresponding to the SM gauge symmetries of U (1), SU (2) 
and SU{3), respectively, qLi, u*j^i, d*j^-, in and e*j^- are respectively the scalar component 
fields of Qli, Uf^^, D'j^-, Ln and £'^j, and $j {j = 1-4) are the scalar component doublet 
fields of the chiral superfields Hj. From W and C^^^^, the Lagrangian is constructed as 



^ =-^^kinctic + -^^gauge-matter " ^- ^^ • IpLj + h.C^ - -{Qaf {'f*T^j<fj , 

+ ^soft ' (2-3) 

where V'Li and ipj represent fermion and scalar component fields of chiral superfields in the 
model. 

There are two Higgs doublets for each quantum number, so that they can mix with each 
other. The Yukawa sector then produces a dangerous FCNC via the scalar boson exchange 
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Table 1. Classification for the charge assignment for the Z2 symmetry in the 4HDSSM. Type C 
and Type D are introduced only when iVj^j are added to the model. 

at the tree level. There are several ways to eliminate such an excessive FCNC. In non- 
SUSY extended Higgs sectors with multi-doublets, a softly-broken discrete Z2 symmetry 
is often imposed [27]. In the general two Higgs doublet model, there are four types of 
Yukawa interactions under such a Z2 symmetry depending on the assignment of the Z2 
charge [29-31]. The other possibility of eliminating the FCNC may be to consider a certain 
of alignment in the Yukawa sector [32], but we do not consider this possibility in this paper. 
In the 4HDSSM, we also impose the Z2 symmetry to eliminate the FCNC. There are two 
types of Yukawa interactions (Type A and Type B) as shown in Table 1, assuming that all 
the Higgs doublet fields receive VEVs. If we introduce additional chiral superfields A'']^j for 
right-handed neutrinos which are singlet under the SM gauge symmetries, possible number 
of the type of Yukawa interaction becomes doubled under the Z2 symmetry, depending 
on the two possible assignment of the Z2 charge for N'^^. We define additional two types 
in Table 1 (Type C and Type D) which correspond to the Z2 odd iV]^^. Under the Z2 
symmetry, some of the Yukawa coupling constants are forbidden for each type of Yukawa 
interaction. For example, in the MSSM-like Yukawa interaction (Type A) 1^ = = Y^' = 
is required, while in the lepton specific one (Type B) we have 5^ = = 1^ = 0. Marshall 
and Sher discussed phenomenology of the Type B Yukawa interaction in the 4HDSSM [28]. 
Notice that the dimensionful parameters are not forbidden as long as the discrete symmetry 
is softly broken [33]. In this paper, we assume that the FCNC is sufficiently suppressed by 
a softly-broken Z2 symmetry. However, we do not specify the type of Yukawa interaction, 
because all the essential results in this paper do not depend on the types of Yukawa 
interaction. The general discussion on phenomenological consequences in flavor physics is 
given elsewhere [35]. 

From the Lagrangian in Eq. (2.3), we can extract the Higgs potential of the model, 
in which each neutral scalar component $i {i = 1-4) receives the VEV. However, because 
Hi and H^ {H2 and H/i) have the same quantum numbers under the SU{2) x U{1) gauge 
symmetries, there are U{2) symmetries in the D-terms [(<?a)^(y'a^Q^¥'/3)^] in the potential. 
By using the U{2) symmetry, we may rotate the fields $1 and $3 and take the basis in 
which only one of the doublets receives the VEV while the other does not. We also can 
perform the same procedure for $2 and $4: see Appendix A for details. Consequently, 
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without loss of generality we can rewrite the Higgs potential as 



^1 



(MDh (Ml) 
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22, 
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(M|)n (M|) 
(M|)2i (M|) 



12 



22 
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(^$^^>2 + $2^$2 - ^1^1 - '^1^1 



+ y {($l$2)(^'^^l) + I^'2)K^^l) + (^;^$2)(^^^'l) + {^'l^'2m'^<^[) 

+ {<i>l<^[){<i>'l<l>,) - {<i>l<^,){<i>'l<^[) + (<I>|c&'2)(<D'j$2) - (<I>|cI>2)(<D'jc&'2)} , 

(2.4) 

where $i (Y = -1/2) and $2 = 1/2) have VEVs, while ^[ {Y = -1/2) and $'2 
(y = 1/2) do not. Throughout this paper, we restrict ourselves in the CP invariant case. 
We thus hereafter neglect all CP violating phases in the dimentionful parameters. 
The rotated Higgs doublet fields $1, ^2 and ^'2 are expressed as 
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-^1 



$5 



V2 
^'2 



, ^'1 



(2.5) 



where the neutral scalar fields can be parameterized as 

V'l = (^^1 +(t>i+ iXi) , V'" = (t^2 + <^2 + iX2) 
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(2.6) 



where the VEVs of these neutral components are given by {(p^) = {1P2) = V2/V^, 

{(fi) = and {ip2) = 0. Introducing 



{V2Gp)-^/'^ ~ 246 GeV 



(2.7) 



and the mixing angle /3, we express vi and u 2 as vi = v cos f3 and V2 = v sin /?. The vacuum 
conditions for the Higgs potential are given by 
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Solving this set of conditions, one can eliminate {Mi)^^, (M|);^^, (M^)]^25 ^-^^d {M2)i2- 

After imposing the vacuum conditions, the mass matrices M^, and for the 

CP-odd, charged and CP-even scalar component states are respectively obtained in the 
basis of ($1, 5 $'ii ^2)- It is however more useful to work the mass matrices of the CP- 
odd scalar bosons and the charged Higgs bosons in the gauge eigenstate basis (the so-called 
Georgi basis) as [36] 



Ml =OlMlOo 
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(2.10) 



Oo 



(2.11) 



with the orthogonal matrix 

/ S/3 0\ 

-S/3 C/3 

10 
\ 1/ 

where we used the abbreviation such as sin 6 = sg and cos 6 = cg. In this basis the massless 
modes, which are Nambu-Goldstone bosons to be absorbed by the longitudinal modes of the 
weak gauge bosons, are separated in the mass matrices. The basis taken here is essentially 
the same as that discussed in Ref. [26]. It is also useful to rotate the mass matrix for the 
CP-even scalar bosons as 
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C2p) 



3 Decoupling property of the extra doublet fields 
3.1 Nondecoupling effects in the large mass regime 

In general, new physics can be tested not only by direct searches but also by indirect 
searches. The indirect searches are performed by precise experiments to find effects of a 
heavy new physics particle on the observables which are well predicted in the low energy 
theory such as the SM. Such new particle effects on the low energy observables usually 
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decouple in the large mass limit after the renormalization calculation is completed. This 
is known as the decoupling theorem proposed by Appelquist and Carazzone [25]. It is 
also known that the decoupling theorem does not hold when the new particles receive their 
masses from the VEV of the Higgs boson. In fact, there is a class of the new physics models 
where nondecoupling effects of heavy particles can appear on the low energy observables. 
For example, chiral fermions such as quarks and charged leptons cannot have the mass 
term because of the chiral symmetry, so that their masses are generated after the chiral 
symmetry is spontaneously broken by the VEV. Therefore, the effect of a heavy chiral 
fermion does not decouple, and it appears as powerlike or logarithmic contributions of the 
mass in the predictions for the low energy observables. Another example is the additional 
scalar fields in extended Higgs sectors whose masses are typically expressed as 



where A' is the coupling constant with the SM-like Higgs boson and M is the invariant 
mass parameter which is irrelevant to the VEV. When X'v"^ > M^, the situation is similar 
to the case of chiral fermions, so that the effect of these scalar bosons do not decouple in 
the large m^p regime. 

The indirect effects of these nondecoupling particles appear in the low energy observ- 
ables at the tree level, or at loop levels such as the electroweak 5", T and U parameters [37] 
and vertex corrections to the SM coupling constants like the 77/i vertex [38] and the hhh 
vertex [39, 40]. In particular, the nondecoupling effects of chiral fermions and additional 
scalar bosons to the triple Higgs boson coupling are known to give quartic powerlike con- 
tributions for the heavy particle masses in the radiative corrections, so that the effect can 
be very significant. As shown in Ref. [40], in the two Higgs doublet model, the deviation 
from the SM prediction can be of order 100 % without contradiction with perturbative 
unitarity [41]. This nondecoupling effect on the triple Higgs boson coupling would be ap- 
plied to realize the strong first-order phase transition [42] which is required for successful 
electroweak baryogenesis [11, 43, 44]. It is obvious that an excessive nondecoupling effect is 
bounded by the theoretical constraints such as perturbative unitarity [45] and triviality [46] , 
because the large mass of new particles directly means a large coupling constant. 

Let us consider the effect of the heavy particles in supersymmetric standard models. 
In general, a SUSY Higgs potential is composed of the D-term, the F-term and the soft- 
breaking term; 



Quartic coupling constants in the potential can come from both the D-term and the F-term. 
In the MSSM, however, because of the multi-doublet structure only D-terms contribute to 
them, which are given by gauge coupling constants. Consequently, the mass of the lightest 
CP-even Higgs boson is determined by the gauge coupling constants and the VEVs at the 
tree level, which is less than mz- A substantial F-term contribution enters into the Higgs 
potential at the one-loop level via the superpotential W ~ —ytUfi^H2QL3 + where yt 
(= V2mt/ f sin/3 ~ (¥11)33) is the top Yukawa coupling constant. The corrected mass of h 




(3.1) 



V = + |Fp + (soft-breaking term). 



(3.2) 
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is described by 

ml ~ ml cos2 2/3 + In (3.3) 

This one-loop correction shows a nondecouphng property in the large mass limit of stops. 
Consequently uih can be above the LEP bound at least when one of the stops is heavy 
enough. There are also contributions from the fx parameter and the soft-breaking At (A^) 
parameter when there is the left-right mixing in the stop (sbottom) sector. Their one- 
loop effects can also be nondecouphng and then can be significant to some extent when 
they are taken to be as large as the scale of the soft-SUSY-breaking mass mgusY- In the 
NMSSM and the MSSM with triplets, nih can be significantly enhanced by the F-term 
contribution [19-21]. Notice that these F-term contributions should vanish in the SUSY 
limit due to the nonrenormalization theorem. These F-term contributions also affect the 
prediction on the other SM observables such as the triple Higgs boson coupling constants 
Xhhh similarly to the case of non-SUSY extended Higgs models which we already discussed 
above [22]. 

On the other hand, a typical example for extended SUSY Higgs sectors without inter- 
actions from the tree-level F-term is that with only multi-doublet structures, such as the 
4HDSSM. In this class of models, if there is no mixing between the light two doublet fields 
and the additional ones, the effects of the extra fields on the MSSM observables become 
suppressed due to the decoupling theorem when the extra doublet fields are heavy, and 
the model behaves like the MSSM. However, nonvanishing effects can appear through the 
B-term mixing when (M|)i2 or (M|)2i grows with {Mf)22 or (M|)22- These effects appear 
at the tree level, so that they would give substantial modifications in the predictions in the 
MSSM for the low energy observables. We stress that these nonvanishing effects due to the 
B-term mixing are not the nondecouphng effects which appear in the large mass limit for 
the new particles when A'f ^ > as a consequence of violation of the decoupling theorem. 
In this sense, we call the nonvanishing B-term mixing effect as the quasi- nondecouphng 
effect. Notice that the quasi- nondecouphng effect only appears in the predicted values in 
the MSSM. It gives modifications in the MSSM predictions such as the masses of h, H 
and and the mixing angle a as well as coupling constants for the MSSM-like Higgs 
bosons. Such an effect, however, does disappear in the predictable SM coupling constants 
of /i77, hWW, hZZ and hhh in the SM-like limit {niA — > oo) according to the decoupling 
theorem. 

Therefore, we would like to address the question of how the extra doublet fields in the 
extended SUSY model can affect the observables which appear in the MSSM, such as the 
mass nifj, {(f) represent /i, H and H^), the mixing angle a, the vertex F^iyy {V = 
and Z; (j)' = h and H) as well as quarks and leptons Y^»ffi {f,f' represent fermions and 
(j)" = h, H, A and H^). Deviations from the renormalized MSSM observable parameters 
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may be expressed as 



MSSM 
2/« „,mMSSM 



~ mT""'" (1 + 



sin2(/3 - a,ff) ~ [sin^(/3 - a)]^^^««^^^ (1 + 5, 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



where (5^, 5s, S^jjiyv oi' ^<t>"ff' represent the effect of the extra heavy scalar fields on each 
observable in the extended SUSY models. In this paper, we study 5h, Sh, 5h± ^iid 5s in 
the 4HDSSM. The MSSM predictions are evaluated at the one-loop level using the on-shell 
renormalization scheme in Ref. [14]. We do not discuss the effect on the Yukawa coupling 
constants in this paper, which will be studied in details elsewhere [35]. 

3.2 Definition of the large mass limit and the decoupling property in the 
4HDSSM 

The soft-breaking mass parameters {M^)ij come from the B-terms in Eq. (2.2). When we 
consider the case with (M|)i2 = (Af|)2i = 0, the mass matrices M^, and are 

block diagonal. The upper 2x2 submatrix in each mass matrix corresponds to that in the 
MSSM; i.e., 2(M|)ii/s2/3 — )• and the other 2x2 submatrix corresponds to that for the 
extra two scalar bosons. They are separated completely in this case. The model effectively 
becomes the MSSM in the large mass limit of the extra scalar bosons. On the other hand, 
in the case with nonzero (.^3)12 or (^^3)21, the masses of the light scalars h, H and 
are modified from the MSSM predictions by the mixing via the B-terms between $1 and 
or between ^2 and ^'2- These effects are expected to be nonvanishing when (M|)i2 or 
(M|) 21 grows with taking a similar value to the 3-3 or 4-4 component in the mass matrices 
such as (M^)22 or (Af|)22. We here discuss these effects in details in the following. 

We start from discussing the CP-odd scalar mass matrix. In order to examine the 
decoupling property of the mass matrices, we further rotate as 



Ml 



/l o\ 
10 

Cg SI 

VO -s-g c-ej 



Ml 



/l \ 
10 
Cq -s 

\0 s-e c-e J 



/O \ 

^^Mllix k'M^ kM'^ 
M2 



S20 
fc'M2 



\0 kM'^ 







(3.8) 



where 



tan 20 



2{M? 



3)22 



{M()22 - (M|)22 + m|c2^ 



(3.9) 



and M , k, k' and r are defined such that 



(M|)i2 



{M'i)2i = K21M 



kuM'^ 



{M()224 + {Mi)22sl + {Mi)22S2e + ^Ca^Cs," 



{Ml)22sl + (M|)22C| - {Ml)22S2-e 



rrir 



-C2/3C2e 



rM^ 



(3.10) 



-10- 



and 

k' = —K21 + —Ki2 , k = -^K21 + —Ki2 • (3.11) 

These parameters are relevant to the extra doublets, then the decoupling limit is taken 
as — )• oo. Here we assume that m\ <C and we treat m\/M'^ as an expansion 
parameter. One of the eigenvalues of M\ should be m\, the mass of the lightest CP- 
odd Higgs boson A, which should coincide with 2{M'^)ii/ S2p in the limit of M — t- oo if 
k = k' = 0. In generic cases, after diagonalizing the mass matrix, {M^)ii is expressed in 
terms of as 

2(M|)„ _ f^Hd + ^-'V + o f 4) U M^^^ . (3.12) 



S20 - I r2 \M^^ , 

The mass eigenvalues for heavier states Ai and A2 are 

~ oiM^ , ~ a2M'^ , (3.13) 

where ai and 02 are given by 



k"^ + + ^(1 + r) - a/ {A;2 + _ _ i)}^ _^ 4/j2^(y, _ 1) 
= 2; 



k"^ + rk''^ + r(l + r) + J{A;2 + rA;'^ - r(r - 1)} + Ak'^rir - 1) 
«2 = ^ • (3.14) 

We note that ai — )• 1 and 02 — )■ r for /c = A;' — 0. 

For the charged Higgs mass matrix, via the similar procedure to the case of the CP-odd 
Higgs bosons, we obtain the deviation in mu± , the mass eigenvalue for the lightest charged 
scalar i?^, from the MSSM prediction as 

mH± = ^(m2^±)MSSM(i + 5^^)^ (3.15) 

where 

1 ml, k^ + k'^r^-C2p{{k^-k'^r^)c^s + 2kk'rs^g} fm\\ 

'^^^ 2m\ + ml, + (1 + A:'2)r2} ^ ^ \AP J ' ^'^^ ^^ 

and {m'^±)^^^^ is the prediction in the MSSM renormalized in the on-shell scheme [14], 
which is simply given by [47] 

(m^,,)™ =m\ + ml. - H^f {m\ + ml,) + TiTA^m\) + U'^Umh), (3.17) 

where (p^) represent the one particle irreducible diagram contributions to the two 
point function of the field (j) at the squared momentum p'^. Masses of the heavier charged 
scalar bosons and Hf^ are obtained as 

± ~ aiM^ , m'i± ~ a2M^ . (3.18) 
^1 ^2 
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Figure 1. The deviation 5h± defined in Eq. (3.15) due to the quasi-nondecoupling effect of the 
B-term mixing parameterized by k and k' in the 4HDSSM. We here take M = 500 GcV, r = 1 and 
^ = 0. The SUSY breaking scale for the MSSM particles is taken to be 1 TeV, and the trilinear 
soft-breaking parameters At and Ai, as well as the /i parameter are taken to be zero. The upper 
figures: 5fj± as a function of for tan/3 = 3 (left) and tan/3 = 10 (right) for k = 1.0, 2.0 and 5.0 
with fixed fc'(= 0.0). The lower figures: 5h± as a function of k for tan/3 = 3 (left) and tan/? = 10 
(right) for k' — 0.0, 2.0 and 5.0 with the fixed niA (= 150 GeV). In all figures, the solid curves 
are the results from the full numerical calculation, while the dotted curves are those by using the 
approximated formula in Eq. (3.16). 



In Fig. 1, we show the numerical results for the deviation 6ij± defined in Eq. (3.15) due 
to the quasi-nondecoupling effect of the B-term mixing parameterized by k and k' in our 
model. The solid curves in the figures represent the results from the full numerical calcu- 
lation, while the dotted curves are those by using the approximated formula in Eq. (3.16). 
The deviation 5ij± turns out to be negative, and amounts to —20 % for a relatively small 
value of ruA- The magnitude of the deviation is smaller for a larger value of niA, but still 
a few times —I % even for = 300 GeV. On the other hand, the deviation is not very 
sensitive to tan /3. We note that the results are insensitive to the details of the MSSM 
parameters such as the soft-breaking mass parameters, the parameter and the trilinear 
Atfi parameters. In fact, when /x and At^i, are varied in the phenomenologically acceptable 
regions, the radiative corrections vary at most from —2 % to +2 %. We have confirmed 
that our results on the one- loop correction in the MSSM agree with those given in Ref. [47]. 
The mass of can be determined with the accuracy of a few percent via the decays of 

—7- Tu and — ?• tb at the LHC [48] , and with the statistical error of less than 1 % via 
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e+e^ —7- H^H~ at the ILC [49]. The mass of A can also be determined with the resolution 
about 2% via the decays A — t- at the LHC, while at the ILC it can be measured 

with the precision 0.2 % via e"'"e~ — t- HA [49]. Therefore, the quasi-nondecoupling effect 
on mH± can be extracted when both mH± and are measured at future collider exper- 
iments. The prediction on mu± (not on 6fj±) in the 4HDSSM is shown in Fig. 5 with the 
comparison of the result in the MSSM. 

Next, the CP-even scalar mass matrix Mfj can also be diagonalized. We first define 
Ml by 

/lOO 0\ /lOO o\ 
10 ,-.,0 10 



M 



H 



Cg 

VO s-s 



M: 



H 



J 



Cq 
VO 



s 

%7 



(3.19) 



and according to the usual mathematical procedure Mfj can be block-diagonalized by 
rotating the basis with an appropriate orthogonal matrix Omh as 



Omh^rOmh 



( 



\ 








where R is defined as 



aiM2 





R 









1 + 4 + A;'2 



+ ' (3-20) 



(3.21) 



The upper 2x2 submatrix coincides to the mass matrix of the two light scalar bosons H 
and h of the MSSM when M — t- oo if /c = A;' = 0. For the case with nonzero k and k\ 
after diagonalizing the 2x2 submatrix by the mixing angle Oeff the mass eigenvalues of 
the CP-even Higgs bosons are obtained as 



2 1 

h 2 



+ ml {clp + R^s%) -M 



l-(l-i?2)si^)} +4m>|.2^i?2 



rrir 



m, 



loop 



m 



H 



nrir 



l-(l-i?2),2^)} +4m>|.i^i?2 

(3.22) 



where AJ^°^ and A'^°'^ represent the one-loop corrections in the MSSM. The masses of 
heavier states H[ and H2 are given by m^, ~ oiM^ |l -|- 0(m^/M^)} and m^, ~ 
02 {1 + 0{m\/M'^)^ . The mixing angle satisfies the relation 



tan(/3 — Oeff) 



•m\c2l3S2l3R 



+ A. 



loop 

tan(/3— a) 



(3.23) 
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Figure 2. The deviation Sh in m/j = mj^^^'^(l + 6h) due to the quasi-nondecoupling effect of 
the B-term mixing parameterized by k and k' in the 4HDSSM, where 

^MSSM ig 

the renormalized 

mass of h. We here take M = 500 GeV, r = 1 and ^ = 0. The SUSY soft-breaking scale for the 
MSSM particles is taken to be ITeV (solid curves) and 2 TeV (dotted curves), and the trilinear 
soft-breaking parameters At and Ah as well as the fj, parameter are taken to be zero. The upper 
figures: 6h as a function of niA for tan/3 = 3 (left) and tan/3 = 10 (right) for k = 1.0, 2.0 and 5.0 
with fixed k'{~ 0.0). The lower figures: 6h as a function of k for tan/3 ~ 3 (left) and tan/3 ~ 10 
(right) for k' = 0.0, 2.0 and 5.0 with the fixed tua (== 150 GeV). 

where ^[^n^^-a) one-loop correction in the MSSM. Notice that ruh and rriH given in 

Eq. (3.22) and tan(/3 — Ocs) in Eq. (3.23) do not depend on the sign of k and k' . We note 
that the effective mixing angle a^s contains information of the B-term quasi-nondecoupling 
effects between $i and or between $2 and ^'2 by k and k', but for m\ <^ the tree 
level formula with the angle a in the MSSM can still hold by replacing a by aoff in a good 
approximation. For example, the coupling constants of the two light CP-even Higgs bosons 
with the weak gauge bosons V {V = and Z^) in the case with nonzero k and k' are 
given by 

^VVh = {c/3(Oh)i2 + sp{Oh)22} = ^ sin(/3 - aeff) (1 + A;^°°p.) , (3.24) 

^VVH = MOh)ii + spiOH)2i} = ^ cos(/3 - aeff) (1 + ^Hvv) ' (3-25) 

where the matrix Oh is given in Eq. (B.5) in Appendix B, and Aj^yy and A^^^y represent 
radiative corrections in the MSSM. Finally, in general, magnitudes of k and k' are not 
necessarily smaller than 1, still it is helpful to deduce the approximate formulae assuming 
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Figure 3. The deviation Sh in mn = m^^^^{l+SH) of the renormalized mass of the second lightest 
CP-even Higgs boson H due to the quasi-nondecoupling effect of the B-term mixing parameterized 
by k and k' in the 4HDSSM. We here take M = 500 GcV, r = 1 and ^ = 0. The SUSY soft-breaking 
scale of the MSSM particles is taken to be 1 TeV (solid curves) and 2 TeV (dotted curves), and 
the trilinear soft-breaking parameters At and At as well as the fi parameter are taken to be zero. 
The upper figures: Sh as a function of ttia for tan/? = 3 (left) and tan/3 = 10 (right) for k = 1.0, 
2.0 and 5.0 with fixed k'{— 0.0). The lower figures: 6h as a function of k for tan/3 = 3 (left) and 
tan/3 = 10 (right) for k' = 0.0, 2.0 and 5.0 with the fixed (= 150 GeV). 



that they are small; 



(3.26) 

'(m^ — m|)2 -|- ArrY^rn^s^p \ 

(3.27) 

tan(/3 - aefr) =[tan(/3 - a)]^SSM I ^ ^ - 2rnl - n^.s^C^ + k'^) ^ 

+o(^)), (3.28) 



where (m|)^^^'^, (rrvjj)^^^^ and [tan(/3 — a)]'^^^'^ are the corresponding parameters eval 



-15- 



0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 







2.0 


tan |3 = 3 
k' = 0.0 


yf^" 


k = 


1.0 






5.0 















100 150 200 250 

mA[GeV] 



300 



0.12 
0.1 
0.08 
0.08 
0.04 
0.02 




-k' 


= 5.0 








= 2.Q/^ 










tan 


3=3 




/k' = 0.0 


TTlA 


= 150GcV 





1 2 


3 


4 




100 150 200 250 

mA[GeV] 



0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




300 









k' = 


5.0'' 




k' = 




tan |3 = 10 




1c' = 0.0 


TTlA = 150GcV 



Figure 4. The deviation 5s in Eq. (3.5) due to the quasi-nondecoupling effect of extra doublet 
fields via the B-term mixing parameterized by k and k' in the 4HDSSM. We here take M ~ 500 
GeV, r = 1 and ^ = 0. The SUSY soft-breaking scale of the MSSM particles is taken to be 1 
TeV (solid curves) and 2 TeV (dotted curves), and the trilinear soft-breaking parameters At and 
Afy as well as the /i parameter are taken to be zero. The upper figures: 5s as a function of tua for 
tan^ = 3 (left) and tan^ = 10 (right) for k = 1.0, 2.0 and 5.0 with fixed k'{^ 0.0). The lower 
figures: 6s as a function of k for tan/3 = 3 (left) and tan/3 = 10 (right) for k' = 0.0, 2.0 and 5.0 
with the fixed tua (= 150 GeV). 

uated at the one-loop level assuming the MSSM. In this paper, we have used the ap- 
proximate one-loop formula given in Ref. [14] in evaluating (m^)'^^^^, {rin?fj)^^^^ and 
[tan(/3 - a)]^^SSM_ 

In Fig. 2, we show the numerical results for the deviation 
6h), where 

^MSSM 

the one-loop corrected mass of h, due to the quasi-nondecoupling 
effect of the B-term mixing parameterized by k and k' in the 4IIDSSM. The SUSY soft- 
breaking scale of the MSSM is taken to be 1 TeV and 2 TeV, and the trilinear soft-breaking 
parameters At and A), and the fx are taken to be zero. It is found that 5h is always positive. 
This is understood from Eq. (3.20). The parameter R is unity for k = k' = 0, and is 
smaller for larger values of k and k' . A smaller value oi R (R < 1) reduces the value of 
the off-diagonal term in Eq. (3.20), which makes the mixing between the first two CP-even 
states weaker. Consequently, the mass difference between h and H becomes smaller than 
the case with the MSSM case with the same value of uia and tan f3. The deviation takes its 
maximal values (6-20 % for tan /3 = 3 and 2-5 % for tan /3 = 10) around the crossing point 
{mA ~ 130-150 GeV) where the role of h and H are exchanged. For larger values of ttia 
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the magnitude of 6h is smaller, but it can be still 3-6 % (about 1 %) at uia = 200 GeV for 
tan/3 = 3 (10). These values are substantial and can be tested by the precise measurement 
of rrih at the LHC (the ILC), where rrih is expected to be determined with about 0.1% [50] 
accuracy at the LHC, while at the ILC it is expected to be measured within less than 70 
MeV [51]) error. The prediction on nih (not on 6h) in the 4HDSSM is shown in Fig. 5 with 
the comparison of the result in the MSSM. We can see that in the 4HDSSM m/j reaches 
its maximal value at a smaller uia than that in the MSSM, although the predicted upper 
bound on the nih is the same in both models. 

In Fig. 3, we show the deviation dn in itt-h = m^^^^(l + Sh), where rrijj^^^ is the 
one-loop corrected mass of H, due to the quasi-nondecoupling effect of the B-term mixing 
parameterized by k and k' in the 4IIDSSM. The SUSY parameters are taken as in the same 
way as Fig. 2. As we discussed, the mixing of the light two CP-even states is weakened 
by non-zero values of k and k', so that mn is smaller than the prediction in the MSSM. 
Therefore, Sh is negative as we expect. The behavior of 6h as a function of ruA and tan (3 
are similar to the case of 5h except for the sign. The magnitude is maximal around the 
crossing point (m^ = 130-150 GeV), and amounts to —18 % (—5 %) for tan/3 = 3 (10). 
At the LHC and the ILC, the mass of H can be determined with the similar precision to 
that of A mentioned in the previous paragraph. The prediction on niH (not on 5h) in the 
4HDSSM is shown in Fig. 5 with the comparison with the result in the MSSM. 

In Fig. 4, we show the numerical results for the deviation 6s defined in Eq. (3.15), in 
which [sin^(/3 — a)]^^^^ is the one-loop corrected mixing factor sin^(/3 — a) evaluated in the 
MSSM. 5s is the net deviation from the MSSM prediction due to the quasi-nondecoupling 
effect of the B-term mixing parameterized by k and k' in the 4HDSSM. The SUSY soft- 
breaking scale of the MSSM is taken to be 1 TeV and 2 TeV, and the trilinear soft-breaking 
parameters At and Ab and the fi parameter are taken to be zero. In the figures, we can 
see that 5s is negative when uia is smaller than the crossing point at uia ~ 130-150 GeV, 
while it is positive for larger niA- The deviation can be as large as 0(10) % (tan (3 = 3) and 
0(20) % (tan/3 = 10) just above the crossing point; i.e., at around niA ~ 140-150 GeV. It 
is rapidly close to unity for larger values of niA- Notice that for larger soft-SUSY-breaking 
scale, a larger 5s is possible. The prediction on sin^(/3 — acff) (not on 6s) in the 4HDSSM 
is shown in Fig. 5 with the comparison with the result in the MSSM. 

4 Conclusions 

We have investigated the decoupling property of extra heavy doublet fields in the 4HDSSM. 
Even without interaction terms from the tree-level F-term contribution such as in the 
NMSSM, significant quasi-nondecoupling effects of extra scalar fields can occur at the tree 
level due to the B-term mixing among the Higgs bosons. We have deduced formulae for 
deviations in the MSSM observables in the decoupling region for the extra heavy fields. The 
possible modifications in the Higgs sector from the MSSM predictions have been studied 
numerically. 

From the results shown in Fig. 1 to Fig. 5, we have found that the quasi-nondecoupling 
effect from the B-term mixing can be significant in the 4HDSSM, which can change the 
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MSSM observables mjj±, ruh, mn and sin^(/3 — a) to a considerable extent. When the 
Higgs boson h is found via the processes of gluon fusion or vector boson fusion at the 
LHC, rrih is expected to be measured very accurately. The correction 5h due to the quasi- 
nondecoupling effect can be much larger than the expected error at the LHC and the ILC. 
Therefore, we conclude that the effect on nih can be measured. The other MSSM Higgs 
bosons H, A and are expected to be discovered at the LHC as long as rriA is not too 
large. The deviations due to the quasi-nondecoupling effect in these quantities may also 
be identified when they are precisely measured at the LHC or the ILC. 

Detecting the deviations from the MSSM predictions on these MSSM observables, the 
MSSM Higgs sector can be tested, and at the same time the possibility of extended SUSY 
Higgs sectors including the 4HDSSM can be explored even when only the MSSM particles 
are discovered in near future at the LHC and at the ILC. A detailed discussion on the 
Yukawa sector is given elsewhere. 
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A Rotation of the basis 

The Higgs potential in our model is given from Lagrangian in Eq. (2.3) by 



+ {Mlh2H^2 + {MlU^^^ + (Ml )*4c|t|,^ + (M2 )44$t|,^ 




+ y {{H^2)iH^l) + iH^.M^l) + (l';5l'2)(«'^*3) + 



(A.l) 



where g and g' are the SU(2)j;, and U(l)y gauge couplings respectively, and 



{Ml)n =(M!)n + |^i2l' + lrt4 

(M!)33=(Af!)33 + |/i32|2 + |^34 



(M2)^3 =(M^)i3 + /ij2/^32 + Ml4^34 > 

(M2)22=(M|)22 + 1^12!' + 1^321' , 
(M2)44 =iMl)^^ + |/ii4|2 + |^34|2 , 

(M|)24 =(M^)24 + /^12/^14 + /^32/^34 • 



(A.2) 
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Because and $3, and $2 ^^'^ ^4 have respectively the same quantum numbers, we 
may rotate the basis by introducing the 2x2 unitary mixing matrices U-^ and U- as 

ft) -ft)-- ft) ■ ft) - ft) ft)- 

By using this degrees of freedom, one may choose the basis where only and have 
VEV's while those of and ^'2 are zero. In this basis, the Higgs potential is expressed by 

^ ^ H ^0 - (1:?;;: la:) ^K^O ^ H - i:t:) (^^^ 
-((*'^o-te:S::::)^-ft)^^-) 

8 V 

+ ^ {{^l^^M^l) + (^l^^)($2^$l) + i^'^^^M^'l) + {^'l<^'2){^'^^[) 

(A.4) 

Hereafter, we reparameterize the parameters as 



^- \{Ml)\, {Ml) J ^' ' \ {Ml)\2 {M^ 



'22 , 



^ /(M^)22 {MD^A _ f (M|)ll (^I)l2 



(M2)*4(M2),J + l(A/|)!2(M|) 



22 . 



/ i?12A^12 ^14/^14^ f^t , (^^3)12 1 (j^ 5^ 

" V^32/«32 ^34/^34/ ^ ^^1)21 (M|)22 ' ' ^ " ^ 

With the above notation, the Higgs potential can be rewritten as in Eq. (2.4). 



B Yukawa Interactions 

In the gauge eigenstate basis, the Yukawa coupling matrices {Yf)ij and {Yj)ij (/ = u,d,e) 
associated with <I>^ and <I>^ (A; = 1 for f = d,e and k = 2 for f = u) are given as 

{Yu)^j ={U+)l,{Yu)^, + {U+)UY:^^J , iYd)^, = iU^)*uiYd)^, + {U^UY^h , 

(y,),, ={U-)uiYe)^j + {UA\2{Yi)^J , 

={U+)*2,{Yu)ij + [UAUYl^h , {Y^h = {U^yM)^J + {U.)l2{Y!i)ij , 
(y,')i,- ={UA*2l{Ye)ij + {uAUYeh ■ (B.i) 

The formulae here are given in the basis where Yu and Y^ are diagonal, i.e., 

(Xu)^, = ^^5^,, ra., = ^^^*, (ye)., = ^^5., , (B.2) 

fsmp fcosp ■' fcosp 
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where m„. , m^. and are the masses of up-type quarks, down-type quarks, and the 
charged leptons, respectively, and Vij is the Cabbibo-Kobayashi-Maskawa matrix. The 
fields fRi and fn (/ = u, d, e) denote the mass eigenstates, and the left-handed component 
of the quark mass eigenstates are embedded into the SU(2) doublets as 



QLi 



'Li 



(B.3) 



From the superpotential given in Eq. (2.1), the Yukawa interactions between the matter 
fields and the neutral Higgs scalar fields are read as follows: 



C 



V2 



— I 



,COt/3(OA)2a - {OA)la) + ^^{OA)4aj URiULjif''^'' 
"^"^^'^ ^tan/3(0^)2a + {OA)la) + ^^^{OHhA dRdLj^l^'' 



V 

I I "^'^''^^-^ (tan/ 



V2 



V 



K0A)2a + (OA)la) + ^^f^ (0^)30 J CRiCLj^''^'''' + h.C. , (B.4) 



where we parameterize the extra down-type Yukawa couplings as = {Y^)ikVjk, CP- 

even and CP-odd Higgs scalar bosons are written as y?^™"^ = {H,h, H[, H2) and ^p'^'^ = 
{z^,A, Ai,A2) with denoting a Nambu-Goldstone mode, and the mixing matrix Oh and 
Oa are defined as 



oImIOh 



(m]j \ 



m^, 

V ' mlj 



/O \ 


m\ 

\0 ' m\j 



(B.5) 



with {OA)la = {OA)al = ha- 

For the Yukawa interactions with a charged scalar fields, they are found as 
C = ( (cot/3(0±)2„ - (0±)i,) + {Yl,\k{0±)Ao\ VkjUmdLjipt 



+ 



;tan/3(0±)2a + (0±)i„) + (yj*)jfc(0±)3a VikULidRjifi 



+ 



^ (tan/5(0±)2a + (0±)i«) + (F^ )j.(0±)3a i^udRj^t + h-c. , (B.6) 
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where the charged scalar bosons are written as (/^^ = {w^ , H^, , with w'^ being a 
Nambu-Goldstone boson, and the mixing matrix 0± is defined as 



/O \ 



m]j^ 
m" 



m^. I 
^2 / 



(B.7) 



with (O±)io = (0±)qi = 5ia. 

In the decouphng region rrij^/M <^ 1, the Yukawa interactions with the MSSM Higgs 
bosons are given as 



C 



(sin(/3 — Oes) + ^cot /3cos(/3 — acs)) 



iyu)^J \k\R 

V2 r 



cos(/3 - acflf) + ^ URiULjh 



(sin(/3 — Qeff) — i? tan /3 cos(/3 — acfr)) 



in: 



V2 



k \Rcos{f3 — acs) + C 



TUe, 6. 



'b-Jl (sin(/3 — aeff) — -Rtan/3 cos(/3 — acs)) 



m , 



^ ,k'\Rcos{f3 - acff) + O ) \ emeLjh 



'cos(/3 — Ocs) — i?cot /3sin(/3 — Oeff)) 

V 

(cos(/3 — agflf) + -Rtan/3 sin(/3 — acfr)) 



+ i^|A:'|fisin(/3 - aes) + O ( ^ ) (iffidiji/ 



^/2 



(cos(/3 — Ocs) + i?tan /3 sin(/3 — acfr)) 



+ 



( 



m. 



^ ,k'\Rsm{^ - acfr) + O j j> eRiCLjE 



— I 



'^Ui^ij T-, . n 

— -^Rcoifi 

V 



{y^\k\R 

V2 



V 



i?tan/3 + 



i2tan/3 + 



V M2 



d-RidLjA 



V2 



k'\R + 0{j^]}eR,eLjA 
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+ |^^^^^^i?tan/3 + (y;)ij|A;'|i?| UueRjH+ + h.c. . (B.8) 

In the SM limit where both M and ruA are enough heavy compared to mz, the contributions 
from Y'j disappear. On the other hand, effects of extra Yukawa contribution can remain 
in the Yukawa couphngs with H, A, and as quasi-nondecouphng effects, even if M 
is much larger than itla- In general, has non-trivial flavor structure and the flavor 
changing processes are enhanced a lot, if there are off-diagonal elements in {Yu)ij^ 0^d)ij 
and {Ye)ij. In order to suppress a dangerous contributions to the flavor changing processes, 
some mechanism to forbid such the off-diagonal elements is necessary unless ruA is very 
large. Discrete symmetries such as Z2 symmetry are often considered. The way of the Z2 
parity assignment are discussed in Sec. II, and the possible types are listed in Table 1. It is 
interesting to study how the extra Yukawa interaction contribute to flavor measurements 
in these types of Yukawa interaction in the 4HDSSM. Further discussions on flavor physics 
are given elsewhere [35]. 
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Figure 5. The values of ,mh,mH and sin^(/3 — aeff) in the 4HDSSM and the MSSM as a 
function of for fc = 1.0, 2.0 and 5.0. The soft-SUSY-breaking scale of the MSSM is set to be 
1 TeV (solid curves) and 2 TeV (dotted curves). The trilinear soft- breaking parameters At and Ai, 
as well as the fi parameter are taken to be zero. The other parameters are takes as M = 500 GeV, 
r = 1, ^ = and k' = 0. Figures in the left column are for tan /3 = 3 and those in the right are for 
tan/3 = 10. 
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